Strong spin-orbit splitting on Bi surfaces 
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Using first-principles calculations and angle-resolved photoemission, we show that the spin-orbit 
interaction leads to a strong splitting of the surface state bands on low-index surfaces of Bi. The 
dispersion of the states and the corresponding Fermi surfaces are profoundly modified in the whole 
surface Brillouin zone. We discuss the implications of these findings with respect to a proposed 
surface charge density wave on Bi(lll) as well as to the surface screening, surface spin-density 
waves, electron (hole) dynamics in surface states, and to possible applications to the spintronics. 
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Recently, spin-orbit coupling (SOC) on surfaces and 
the resulting splitting of the surface-state bands has 
attracted considerable attention. While it is a well- 
established fact that the reduction of coordination at sur- 
faces and in thin films can lead to pronounced magnetic 
effects, the discovery of a small splitting in the band of 
the sp surface state on the non-magnetic Au(lll) surface 
and its interpretation as being due to SOC by LaShell 
et al. came as a surprise. More sophisticated angle- 
resolved photoemission (ARPES) investigations and cal- 
culations have meanwhile confirmed this splitting Q and 
the combination of the experimental results with first- 
principles calculations do indeed proof that the SOC is 
causing it 0, Later, larger SOC-induced splittings 
were found on the Li-covered surfaces of W and Mo 
and the predicted difference in spin-orientations for 
H on W(110) was confirmed experimentally using spin- 
resolved ARPES 5] . Since these surface states contribute 
only very little to the density of states at the Fermi level, 
the observed spin-orbit or Rashba splitting of these states 
will not show up in transport phenomena. On the other 
hand, surface states of a semimetal would give a promi- 
nent contribution |6|, L| which could make these systems 
interesting for applications in the field of spintronics. The 
surfaces of the semimetal Bi seem to be ideal to advance 
our understanding of SOC on surfaces and how it mani- 
fests itself in experiments. Of particular interest are the 
influence of the SOC on the electron-phonon coupling 
electron and hole dynamics 0, and the possible for- 
mation of surface charge (spin) density waves. The oc- 
currence of strong SOC in low-dimensional structures of 
non-magnetic materials could also have applications like 
spin-filter devices. 

ARPES measurements of the Fermi surface (FS) and 



surface states were recently performed by Ast and Hochst 
for Bi(lll) H3. They interpreted the obtained FS in 
terms of two different surface bands which are not de- 
generate at the f point. Based on this electronic struc- 
ture they proposed a possible mechanism for the forma- 
tion of surface charge density waves (CDW) on Bi(lll) 
Agergaard et al. Q measured surface states and 
the FS on Bi(llO). They pointed out that these surface 
states should be completely non-degenerate because of 
spin-orbit splitting but, as we show below, the splitting 
is so large that an easy identification of the spin-split 
bands was not possible. 

Bulk Bi is a semimetal where the strong spin-orbit in- 
teraction essentialy accounts for the existence of the hole 
FS at the T point 12] but docs not lead to any lifting 
of the spin degeneracy because of inversion symmetry. 
However, the high atomic number of Bi (83) and the 
pronounced splitting in the atomic p levels (the atomic 
P3/2-P1/2 splitting in Bi ( 1.5 eV) is three times stronger 
than in Au (0.47 eV) [ljj) should lead to an observable 
splitting of the surface state bands. 

In this Letter we show that indeed the surface states 
on low-index surfaces of Bi exhibit a spin-orbit splitting 
of the bands which is by far stronger than any case re- 
ported so far. We prove this by combining the results of 
first-principles calculations with high-resolution measure- 
ments of the electronic structure by ARPES. The results 
of the calculation agree well with experiment but only if 
the SOC and hence the removal of the spin degeneracy is 
taken into account. We find that the SOC induced split- 
ting is an essential ingredient for the description of the 
electronic structure: it profoundly changes the surface- 
state dispersion and the corresponding Fermi surfaces on 
all the Bi surfaces of interest. In particular, it is respon- 
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sible for the existence of the six FS hole lobes in the f M 
symmetry directions for Bi(lll), and it also leads to cor- 
rect dimensions of the electron FS hexagon around the f 
point. Our relativistic calculation demonstrates that de- 
spite the existence of nesting at the electron FS hexagon 
of Bi(lll) the formation of a CDW [HI appears to be 
improbable since this nesting couples the states with the 
same energy but different spin, e(k, f) = e(— k, 1). It 
could lead rather to the formation of a spin-density wave 
and not to a CDW. 

Here we present results for three surfaces: Bi(lll), 
Bi(llO) and Bi(100). The calculations have been per- 
formed by using the full-potential linearized augmented 
planewave method in film-geometry as implemented 
in the FLEUR program and local density approxima- 
tion for the description of exchange-correlation poten- 
tial. SOC is included self-consistently as described in 
Ref. [13] • All the Bi surfaces were simulated by a 22 
layer film embedded in vacuum. One side of the film was 
terminated with hydrogen to avoid interaction between 
the surface states of the two surfaces of the film. The 
H atoms were placed in a distance of » 2 A from the 
Bi. On the other side of the film, the termination was 
chosen such that the interlayer distance between surface 
and subsurface atoms was the shorter one of two possible 
terminations. For the calculations a planewave cutoff of 
-f^max = 3.4 (a.u.) -1 was used and the surface Brillouin 
zone (SBZ) was sampled with up to 121 k||-points. 

We also show surface states measured with angle- 
resolved photoemission spectroscopy. The experiments 
were performed at the SGM-3 beamline of the syn- 
chrotron radiation source ASTRID in Aarhus 0|. All 
surfaces were prepared from mechanically polished single- 
crystal surfaces which were cleaned in situ by cycles of 
Ne ion bombardment and annealing to about 473 K. This 
resulted in well-ordered and clean surfaces as judged by 
low energy electron diffraction and Auger electron spec- 
troscopy, respectively. The total energy resolution for 
the data shown below is better than 35 meV. The angu- 
lar resolution of the analyser is ±0.5°. The samples were 
cooled to approximately 30 K. 

Before discussing our results, we briefly explain the 
symmetries relevant for the spin-orbit splitting of elec- 
tronic bands. Time reversal symmetry requires that 
e(k, 1) — e(— k, !). This means that if one has a sur- 
face state at k with a binding energy e and a spin j 
then there must also be a state at — k with the same 
energy but spin J,. This has to be combined with the 
usual space group symmetry. The combination has two 
consequences. First, if the space group contains inver- 
sion symmetry (e(k, |) = e(— k, j)) the bands are obvi- 
ously doubly degenerate. This can happen in the bulk 
but not at the surface. Second, in the case of surface 
states, the splitting has to be zero at some special points 
of the surface Brillouin zone. The application of time- 
reversal symmetry alone yields that one of these points 
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FIG. 1: Surface states of Bi(lll) calculated without (black) 
and with (red) spin-orbit splitting included. The shaded areas 
show the projection of the bulk bands obtained without (vio- 
let) and with (yellow) SOC and their superposition (brown). 



is the centre of the zone T for which k|| = 0. This is 
indeed observed in the dispersion of the spin-orbit split 
states on Au(lll) 0. The combination of time-reversal 
symmetry with translational symmetry dictates that the 
spitting must also be zero for any point which is situated 
half-way between two T points. 

Figdshows the electronic structure of Bi(lll) together 
with the projected bulk band structure for the (111) sur- 
face calculated with and without SOC. The projection 
was calculated using the tight-binding model of Ref. ^3 • 
For clarity, we only show the surface states which are lo- 
cated on the clean surface of the slab. In the case without 
SOC, we find a parabolic f surface state located in the 
nonrelativistic energy gap. Around T this surface state 
band gives an electron FS hexagon. When the SOC term 
is included it results in a spin-splitting of the surface state 
in all the symmetry directions and leaves it degenerate 
only at T and at M. The latter is expected because M is 
a high-symmetry point on the SBZ boundary which lies, 
in contrast to e.g. K, in the middle between two T points. 
Around F this relativistic surface state is degenerate with 
bulk states and shows less clear surface character. The 
lift of the spin degeneracy leads to radical change of the 
surface FS: 1) The radius of the FS hexagon is smaller by 
30% compared to the nonrelativistic calculation; 2) In the 
TM symmetry directions the hole lobes are formed. An- 
other remarkable feature of the Bi(lll) surface electronic 
structure is the very strong anisotropy of the spin-orbit 
splitting: it is «0.2 eV in the TM direction and even 
more in the FK one. 

The occurence of strong spin-orbit splitting is con- 
firmed when the calculations are compared to the exper- 
imental results. Fig. shows the calculated electronic 
structure for the three surfaces together with experimen- 
tal data. In Fig. [2a) we compare the results for Bi(lll) 



3 



and find excellent agreement for the two split surface 
states near T. The experimental results also agree with 
recently published data for Bi(lll) ^1 but the two spli t 
bands appear better resolved here. Here, as in Ref. 101, 
the intensity of both surface states strongly decreases 
close to f . This is most likely due to the overlap with 
the projected bulk band structure. The surface states 
are no longer genuine surface states but surface reso- 
nances which penetrate much more deeply into the crys- 
tal and give a lower photoemission intensity. Although 
it is therefore not simple to decide if the states are in 
fact degenerate at T or not, we find no evidence for the 
latter. This discrepancy with the data of Refs. 0, 0] 
is most likely due to a sample misalignment which later 
was found to have been present 11]. Our interpretation 
is in disagreement with the result from Ast and Hochst 
who used a Bi(lll) bilayer to simulate the surface state 
dispersion of Bi(lll) with a tight-binding model 0|. In 
their case, the calculated bands match the experimental 
dispersion beyond the T point when the bilayer thick- 
ness is increased by 70% with respect to the bulk value 
and when the SOC strength was reduced to 13% of the 
experimental value 1.5 eV 13]. The bilayer, however, 
always has inversion symmetry and therefore this calcu- 
lation yields two spin-degenerate bands near the Fermi 
energy that do not cross at T. 

Fig. Efb) shows the situation near f of Bi(110). The 
experimental results have already been published else- 
where 0. In the theory one can clearly see that, as in 
Bi(lll), the surface state on Bi(110) is degenerate at 
r and splits into two surface states along the symmetry 
lines with one electron per k-point. In contrast to Bi(lll) 
this surface state is unoccupied at f and has negative ef- 
fective electron masses that lead to the formation of the 
hole FS pocket around f 0. This specific behavior of 
the surface state bands is also responsible for the forma- 
tion of the electron FS pocket between Xi and Mi and 
the hole pocket at Mi |?J. In the experiment only the 
lower branch of the spin-orbit split state can be observed 
as it enters the occupied states. Such a situation can be 
highly confusing because the band could be mistaken for 
a simple parabolic hole pocket. 

The scenario of a very steep band and a flatter one near 
r can also be found on Bi(100). This is shown in Ref. [2(1 ] 
and is therefore not presented here. Instead, Fig. [^c) 
shows the situation near the M' point. This point is the 
M point of the quasi-hexagonal SBZ of Bi(100) which is 
not lying on the mirror plane of the SBZ (see Ref. [HI)- 
Since all M points of a (quasi) hexagonal SBZ fulfill the 
criterion of lying exactly in the middle of the line joining 
two f points, we also expect a degenerate surface state 
here. As Fig. [21(c) shows, this is indeed the case. In fact, 
here the bands close to the high symmetry point are so 
steep that the dispersion can not be resolved in the ex- 
periment. In the rest of the SBZ the agreement between 
experiment and calculation is more difficult to find. This 



is due to the deep penetration of the surface states into 
the bulk, such that even a calculation with a 22 layer film 
cannot completely avoid the interaction between the two 
surfaces. Details of the electronic structure of Bi(100) 
will be published elsewhere |2l) . 

The spin-orbit splitting obtained for the Bi surfaces 
is a few times bigger than that of the surface state on 
Au(lll) 0,13], which is to be expected since the atomic 
spin-orbit splitting in Bi is three times larger than that in 
Au. However, the character of the spin-orbit splitting in 
the present case is different from that in Au(lll). For the 
latter, the electronic structure can be described by a free 
electron-like two-dimensional surface state at T. In such a 
case the spin-orbit interaction can be treated by adding a 
so-called Rashba term to the non-relativistic Hamiltonian 
[22| . This leads to a splitting of the surface state which 
is linear in k. In the case of Bi the surface states are not 
free electron-like and they are distributed over the whole 
SBZ, therefore the spin-orbit splitting shows a much more 
complex behavior which can only be revealed by first- 
principle band structure methods. 

The SOC-induced splitting should have some impor- 
tant consequences for the physical properties of the Bi 
surfaces, in particular, for the screening. In the Lindhard 
picture of screening, the susceptibility x(q) is essentially 
given by an integral over all processes where an elec- 
tron hops between an occupied state and an unoccupied 
state separated by q. In a two-dimensional situation this 
type of screening can lead to a CDW-type instability only 
when there are 'nested' elements of the FS, separated by 
q. Such a situation exists for Bi(lll) where the FS of 
the electron pocket around f is hexagonal 0] . Ast and 
Hochst have recently shown 0] that the leading edge of 
the energy distribution curves at the Fermi level cross- 
ing shifts discontinuously as a function of temperature, 
indicating the formation of a CDW. However, when we 
take into account the spin of the states involved in the al- 
leged formation of the CDW, the electron hopping across 
the FS would have to undergo a spin-flip because of the 
split nature of the bands. This makes the occurrence of 
a CDW very unlikely. 

The spin-orbit splitting in surface bands on the Bi sur- 
faces can also have drastic consequences for electron and 
hole dynamics in surface states. In particular, the sur- 
face response function should include all the spin-flip pro- 
cesses between the split surface bands with different spin 
direction. It can lead to the formation of surface spin- 
density waves even in cases when the nesting at the sur- 
face FS does not occur. The spin-orbit splitting should 
also lead to different hole (electron) lifetimes in surface 
states compared to that for the non-split surface state. 
This is due to both the surface response function that 
now includes s pin- flip processes and to a different phase 
space factor |fj.l23j . The surface state spin-orbit splitting 
can also affect the electron-phonon (e-ph) coupling on 
the Bi surfaces. Strong e-ph coupling in surface states on 
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FIG. 2: Calculated and measured electronic structure in the vicinity of two high symmetry points on three surfaces of Bi. (a) 
F on Bi(lll), (b) r on Bi(llO), and (c) Mi on BUlOO). The small black dots are the projected bulk band structure calculated 
using the tight-binding model of Liu and Allen |T^|. The red filled circles are the calculated surface state energies, thin red line 
is a guide to the eye. The photoemission intensity is linearly scaled from dark blue (minimum) to white (maximum). 



Bi(100) was already discussed by Gayone et al. |20j who 
assumed direct interaction between electron and phonon 
systems. Here we would like to note that the existence 
of spin split surface states also permits a spin-wave me- 
diated e-ph interaction. 

The lack of inversion symmetry and the large spin- 
orbit splitting holds also at Bi-insulator interfaces where 
spin-orbit split surface states turn to interface states. We 
speculate that ultra-thin films of Bi covered with an in- 
sulator become much more effective spin-filters and spin- 
manipulators for injected spin-polarized electrons than 
semiconductors explored in spintronics because the rel- 
evant fc-vectors and the difference of the fc-vectors for 
spin-up and spin-down electrons are much larger than 
for semiconductors. 

In conclusion we have presented the first ab initio cal- 
culation of the relativistic surface electronic structure for 
Bi(lll), Bi(llO), and Bi(100). We have shown that the 
SOC term leads to a strong and anisotropic splitting of 
the surface state bands that profoundly modifies the dis- 
persion of the surface states and the surface FS. The 
calculated results are in good agreement with the ex- 
perimental data if the spin-orbit splitting is taken into 
account. We have discussed the possible effect of the 
SOC interaction on the surface response function and 
new mechanisms of electron and hole decay in the Bi 
surface states in terms of electron-electron and electron- 
phonon interactions. We also discussed the possible use 
of this effect for spintronic applications. 
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